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Abstract 

We perform high quality, first principles calculations of the prop- 
erties of Pb and Tl isolated monolayers. Among these, we consider 
the equilibrium lattice constant, the two dimensional compressibili- 
ties and the electronic density. Comparison is made with previous 
results obtained using more simplified models. The present results 
represent an improvement concerning the calculated compressibilities; 
these remaining still lower than the measured values. We speculate 
that the latter could be due to some corrugation of the monolayer, not 
considered in the present modeling. 

Keywords: work function, local pseudopotentials, density func- 
tional calculations. 



* Corresponding author. Fax: +54 51 334174; E-mail: eleiva@fcq.uncor.edu 



1 



1 Introduction 



The application of in-situ surface X-ray scattering techniques to the study of 
electrochemical systems has contributed to the elucidation of the structure of 
metal monolayers adsorbed at underpotential Jl], ||, [3|, |], [5|, ||. Among these 
systems is particularly interesting the case of sp metals adsorbed at under- 
potential on dense substrate faces like Ag(lll) and Au(lll). In these cases, 
the adlayer presents an incommensurable structure, whose lattice constant 
varies with the applied electrode potential. The monolayer is built with a 
lattice constant which is smaller than that expected from the bulk structure. 
This contraction of the adsorbed monolayer has been theoretically explained 
by Leiva and Schmickler [12| by means of a jellium model with pseudopo- 
tentials, using the density functional formalism to write down the energy of 
the system. In the case of the detachment of a single slab from the bulk, it 
was found that the electronic density of the isolated layer expanded in the 
direction perpendicular to the lattice plane. This relaxation, concomitant 
with a decrease of the nearest neighbor distance between the ion cores in the 
slab, occurs because the absence of the ions of neighboring planes sitting on 
hollow sites reduces the electronic pressure in the direction perpendicular to 
the slab. In the case of the experiments, it was found that when the potential 
is shifted cathodically, more atoms are packed into the adsorbed monolayer, 
with the concomitant decrease of the lattice constant. From this experimen- 
tal information, the 2D isothermal compressibility Kid is estimated to be 
in the range of 1.2 to 1.9 A 2 /eV. A first estimation of this compressibility 
in terms of a free electron model yielded values in the order of 0.5 A 2 /eV, 
that is, considerably lower than in the experiment. Recent calculations by 
Leiva and Schmickler [|nj using a more sophisticated jellium-pseudopotential 
model also rendered values which were low as compared with the experiment. 
It is the purpose of this work to perform high quality, first principles calcu- 
lations for some of the systems addressed above, which should clarify if the 
relatively low theoretical compressibilities are due to simplifications in the 
model or to some other aspects so far not taken into account. We consider 
here the cases of Pb and Tl, which are the systems for which a wealth of 
information exists. 
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Figure 1: Schematic illustration of the supercell geometry employed in this 
work in order to represent a single layer of a metal isolated in vacuum. 

2 Model and calculation method 



In these preliminary calculations we only consider a metal slab isolated in 
vacuum, which we contract isotropically, thus obtaining the energy as a func- 
tion of the distance between nearest neighbors. 

We performed two sets of calculations. One of them was performed by 
means of a computer code developed by one of us |H |, which solves an effec- 



tive one dimensional Schrodinger equation in the presence of averaged local 
pseudopotentials |L5| within the density functional approach [[18] . Within this 
model, which will be referred as 1-D model onwards, the electronic density 
is a function of the distance perpendicular to the metal surface. The other 
set of calculations was performed with the complex computer code fhi96md 
f|, developed at the Fritz-Haber-Institut by Schemer and coworkers. In this 
case, the geometry actually corresponds to that of a periodic supercell as 
schematically shown in Figure 1 and the pseudo electronic density contains 
the full three dimensional features outside a selected core region. We shall 
refer to the results obtained by this method as those of the 3-D model. The 
geometry of the nuclei is contained in a supercell, which is periodically re- 
peated as a lattice. Thus, the coordinates of a nucleus or its periodic image 
FL will be 
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Ri = n + R 



where Tj represents the position of nucleus "z" inside the supercell and R is a 
lattice vector. The compact surface structure subject of the present studies 
was represented by a simple hexagonal lattice. In order to ensure that the 
results of the calculation accurately represent an isolated surface, the vac- 
uum region was varied and the convergence of the energy was monitored. We 
found that a distance of five times the distance between nearest neighbors 
is adequate for the present purposes. Since the features of the computa- 
tional method are extensively described in the reference given above |J, we 
just give here a short comment on those aspects that are relevant for the 
present application. This program is also based on the density-functional 
formalism [pl|, where the variational problem of a many-particle Schrodinger 
equation is transformed into a variational problem of an energy functional 
E[n(r)}, where we represent the electron density of the system with n(r), in 
short n. The energy equation of the electronic system is then usually written 
as: 



nuc—nuc 



E[n] = T s [n] + E H [n] + E e ~ nuc [n] + E xc [n] + E 

where T s [n] is the functional describing the kinetic energy of a system of 
non interacting electrons with density n(r) , E H [n] is the Hartree energy, 
calculated from the corresponding potential V#(r) and E xc [n] is the so called 
exchange and correlation energy. In the present case we have used the so 
called local density approximation (LDA) to the exchange and correlation 
energy, which applies locally the results for the homogeneous electron gas 
obtained by Ceperley and Alder || in the parameterization of Perdew and 
Zunger |10[. The remaining terms correspond to the electron-nuclei (E e ~ n ) 



and nuclei-nuclei (^E nuc ~ nuc ) electrostatic interaction. The electronic density 
n{r) is obtained through the self-consistent solution of the corresponding 
Kohn-Sham equations: 

-iv 2 + V ext (r) + V h (t) + K. c (r)] ^(r) = e^(r) 

where V ext , Vh and V xc are the external potential, the Hartree and the exchange- 
correlation potentials respectively, which are given by: 



Vff(r) = J dv'- 



Till" 
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and 



V xc r = T-rk 



Thus, n(r) is given by: 

occ 

n(r) = l^(r)| 2 

where /(z) is the occupation number of state i. Instead of considering the 
whole electronic system for the calculation of the properties of the slab, 
the effect of the core electrons and the nuclei on the valence electrons were 
replaced by suitable pseudopotentials, constructed according to the schemes 
of Hamann |16] , which were generated and tested for their transferability (17 



using the code developed by Fuchs et al . |B1 This involves checking for the 
absence of ghost states, which may appear as consequence of using Kleinman 
and Bylander |2D[ fully separable ab initio potentials and monitoring the 
logarithmic derivatives of the solution radial Schrodinger equation at energies 
close to that of the selected reference state. 



3 Chemical potential. Equilibrium condition. 
Compressibility 

When referred to the bulk deposition potential of metal A, the potential 
$(©) at which a certain coverage of A on a substrate S is obtained can be 



written as 21 



*(Q)=IJIA,A-HA,S[ Q ] C 1 ) 

where ^a,a an d /-<vi,s[@] denote the chemical potentials of the A atoms in 
the bulk phase and when are adsorbed on S respectively. We stress the fact 
that this latter quantity is a function of the coverage degree. In the usual 
electrochemical potentiostatic experiment, a given potential is usually set 
and a coverage degree is thus obtained. The calculations follow the reverse 
order. In this certain configuration (coverage degree, nearest neighbor 

distance) of the system is considered, and the chemical potential //^sfO] (and 
thus the potential) can be calculated. 

For the discussion below, we remind the fact that we are calculating the 
ground state properties of an electron gas (T=0), so that our free energies 
are obtained from energy calculations of the system. If we denote with Fa,s 
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the free energy of the substrate/adsorbate system, /ia,s[@] can be calculated 
from: 

^[8] = (2) 

where Na represents the number of adsorbed atoms of type A and constant 
temperature, area and volume are assumed for taking this derivative. Eq. § 
can be also written in terms of the free energy per unit surface F A s = Fa,s/ S 
and the atomic surface a: 

dF A , 3 /S dF\ s 2 dFl s 
^ S[Q] = -dNjS =WJa) = ~ a ^r (3) 

A simple relationship between the chemical potential /iA,s[@] arid the binding 
energy per atom E S A s can thus be obtained by considering that in the present 
case F AS = E S A s . Eq. [3] turns into: 

haj[Q] = -« 2 ^f ( 4 ) 

This equation could also have been formulated in terms of the binding 
energy per adsorbate atom E bind , which can be calculated from E AS and the 
corresponding energy per unit surface of the pure substrate Eg according to 

E Hnd = a (EXs ~ E s s ) (5) 

In the present calculations E S A s will be considered to be equal to the 
energy per unit surface of the isolated slab, since the contribution of the 
substrate will be ignored. 

Using now eqs. f| and we obtain : 

AU,s[@] = E bind - a dE g™ d (6) 

All chemical potentials and binding energies reported in this paper have 
been referred to the energy of the ion cores and the valence electrons infinitely 
separated from each other. The corresponding values referred to the energy 
of neutral atoms can be obtained by adding the proper ionization energies. 

We now turn to analyze the equilibrium condition in order calculate the 
lattice constant of the adsorbed monolayer. The pressure of a thermodynamic 
system can be calculated from: 
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where E and V are the energy and the volume of the system respectively. 
The constant entropy restriction was left out because we are considering the 
zero temperature case. In the case of three-dimensional solids, setting in 
this equation P = provides a way of predicting the density of a metal if an 
expression for E is available [p2 ]. In the present case, where we are considering 



the isotropic compression of an adsorbed monolayer, it is suitable to state the 
problem in terms of the the extensive coordinate and its conjugate 

force, usually called spreading pressure [p3fl , that we denote with <fi. Thus, in 
the zero temperature case the analog of eq. ^ becomes: 

<~i§ 

where constant temperature and volume are assumed for taking this deriva- 
tive. Thus, we shall employ the condition: 

rn^d = (9) 

to predict the equilibrium lattice parameter of the monolayer. 

The 2D isothermal compressibility of the Tl monolayer was calculated 



according to the usual definition [23 




*2D = ~ hH ( 10 ) 



which is the same as that employed by Toney et al. []7| to get the experimental 
values from the dependence of a nn on the applied potential. 



4 Results and discussion 

In addition to the logarithmic derivatives and the hardness monitoring, a 
further test of the transferability of the pseudopotentials that is usually per- 
formed consists in calculating the equilibrium lattice constant and bulk mod- 
ulus of the metal considered. We performed these calculations for Pb and Tl, 
obtaining the results given in Table 1. As can be seen from the comparison 
with the experimental values found in the table the agreement is very good, 
not only in the lattice constant but also in the bulk modulus. 

Energy calculations as a function of the nearest neighbor distance a nn 
for Pb and Tl isolated slabs are shown in Figure 2. The a nn values at the 
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Figure 2: Energy as a function of the distance between nearest neighbors for 
single Pb and Tl monolayers isolated in vacuum. Results obtained with local 
pseudopotentials, 1-D model a)Pb, b)Tl. Results obtained with non-local 
pseudopotentials, 3-D model, c)Pb, d)Tl. 
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Figure 3: Average electronic density, plotted as a function of the distance 
in the direction perpendicular to the plane of the metal slab for (a) Pb and 
(b) Tl. Broken line: 1-D model, full line: 3-D model. 

minimum, along with the compressibility K 2 d obtained through numerical 
differentiation of these results are shown in Table 2. For the sake of compari- 
son, we show in both cases the results obtained by means of the self-consistent 
calculation employing local pseudopotentials and the one-dimensional solu- 
tion of the effective Schrodinger equation. The curves obtained with the local 
pseudopotentials show a steeper increase of the energy in the neighboring of 
the energy minimum, which is the reason for the outcoming relatively low 
compressibilities. As discussed below, the stiffness of the one-dimensional 
model may be understood in terms of a lack of relaxation of the electronic 
density in the direction parallel to the surface. In this respect, the introduc- 
tion of the real three dimensional structure of the electronic density brings a 
considerable improvement over the calculated compressibilities. 

Figure 3 illustrates some differences concerning the electronic density that 
appear between the 1-D model and the more sophisticated 3-D model calcula- 
tions. In the case of the 3-D model, the electronic density was averaged in the 
direction parallel to the surface (say, the x-y plane) and a one-dimensional 
electronic density was obtained according to: 



It is observed that while the nonlocal pseudopotentials concentrate more 
electronic density on the center of the metal slab, the Ashcroft -empty core- 




I Jn(x,y,z)dxdy 



(11) 
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(c) 

Figure 4: a) Contour lines of the electronic density within a Pb surface 
unit cell, b) Contour plots of the electronic density in a plane perpendicular 
to the surface unit cell, containing two nearest neighbor atoms, c) Plot of 
the electronic density in the line joining two atoms in surface unit cell as 
shown in the inset p mean indicates the average electronic density of bulk Pb. 
Thick point dashed lines in the contour plots show the isoelectronic line 
corresponding to the value of the average valence electronic density of the 
bulk metal. 
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Figure 5: Nearest neighbor distance between the atoms in the monolayer as 
a function of the chemical potential, a) Tl, b) Pb. 



pseudopotentials show a depletion there. On this basis only, and thinking 
in terms of a free electron gas, it could appear that the Ashcroft pseudopo- 
tentials should yield a higher compressibility, since the electron gas with the 



lower electronic density should yield the higher compressibility |24]. How- 
ever, the results of the computations can be understood if we consider the 
electronic density in the x-y plane, as shown in Figure 4a, where we have 
plotted the contour lines of the electronic density within the layer. Figure 4b 
shows similar contour plots for a plane perpendicular to the x-y plane, con- 
taining two nearest neighbor atoms. Figure 4c shows a plot of the electronic 
density in the line joining two atoms in surface unit cell as shown in the 
inset. It can be appreciated that the electronic density in-between the atoms 
is lower than the average electronic density n, thus producing a concomitant 
higher compressibility of the layer as compared with that of a homogeneous 
electron gas with density n. 

Figure 5 shows the nearest neighbor distance as a function of the chem- 
ical potential, calculated according to the formulation of the previous sec- 
tion. The chemical potentials are referred to the chemical potential of the 
layer for = 0. These curves are the theoretical analog of the experimen- 
tal nearest-neighbors vs potential curves. Note that the Tl curves show a 
larger curvature, a fact which is in qualitative agreement with the experi- 
mental findings. These curves should actually not be attainable from the 
electrochemical experiment for (/i — ji eq ) > 0, since this would correspond to 
a negative (stress) of the monolayer, which would produce an exponential 
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increase of the number of defects and thus the breakdown of the monolayer 
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5 Conclusion 

While the present theoretical compressibility values are in qualitative good 
agreement with other previous estimations and show an important improve- 
ment towards agreement with the experiment, they are still too low as com- 
pared with the measured values. Although the presence of a substrate may 
deliver some contribution to the compressibilities, we expect that an adsorbed 
monolayer should be even less compressible than an isolated one, because of 
the contribution of the valence electrons of the substrate. Thus, we think 
that the relatively high compressibility values measured may be related to 
some other aspect of the experiment so far not taken into account. A small 
corrugation of the monolayer under compression could explain this fact. This 
should be considered in the future modeling of these systems. 
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Table 1: Comparison between the theoretical predictions and the exper- 
imental results for the lattice constant and the compressibility of the bulk 
metal. The experimental results were taken from ref. [p5[| . We find a theo- 
retical c/a ratio for Tl of 1.6; the experimental value is 1.599. 1.599. 



Metal 


&(th)/A 


a(exp)/A 


B(th)xl0- 12 dyn/cm 2 


B(exp)/xl0 12 dyn/cm 2 


Pb 


4.98 


4.95 


0.530 


0.488 


Tl 


3.47 


3.46 


0.344 


0.359 
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Table 2: Nearest neighbor distance values (a™ n ) and 2D compressibilities 
obtained at the minimum of the energy vs nearest-neighbor curve. The 
theoretical values correspond to the one-dimensional (thl) and to the three- 
dimensional (th2) calculation respectively. The experimental values (exp) 
were taken from ref. || and correspond to the adsorption of the metals on 
Ag(lll) and Au(lll) surfaces. 



Metal 




ar(^2)/A 


a nn (exp)/A 


Pb 


3.314 


3.271 


3.33-3.43 


Tl 


3.349 


3.212 


3.39-3.50 



Metal 


K%?(thl) 




«2d(exp) 


Pb 


0.48 


0.53 


1.52-1.88 


Tl 


0.28 


0.72 


1.25-1.69 
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